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ABSTRACT 

We investigate whether a single population of first stars could have influenced both the metal en- 
richment and rcionization of the high-rcdshift intergalactic medium (IGM) , by calculating the generated 
ionizing radiation per unit metal yield as a function of the metallicity of stellar populations. We examine 
the relation between the ionizing radiation and carbon created by the first stars, since the evidence for 
the widespread enrichment of the IGM at redshifts z ~ 3 — 4 comes from the detection of C IV absorption. 
We find that the number of ionizing photons per baryon generated in association with the detected IGM 
metallicity may well exceed that required for a late hydrogen reionization at z ~ 6, by up to a factor 
of 10-20 for metal- free stars in a present-day initial mass function (IMF). This would be in agreement 
with similar indications from recent observations of the microwave background and the high-z IGM. In 
addition, the contribution from intermediate-mass stars to the total metal yield, neglected in past works, 
substantially impacts such calculations. Lastly, a top-heavy IMF is not necessarily preferred as a more 
efficient high-z source of ionizing radiation, based on nucleosynthetic arguments in association with a 
given level of IGM enrichment. 



1. INTRODUCTION 

The first stars are expected to be cosmologically signif- 
icant from at least two sets of observations of the high- 
rcdshift (z) intergalactic medium (IGM): (1) hydrogen 
reionization at z > 6 (Becker et al. 2001), and (2) the 
persistent metallicity (Z ~ 0.003Z©) seen in Lya clouds 
up to z ~ 5 (Songaila & Cowie 1996; Songaila 2001). In 
principle, there should be a strong correlation between the 
integrated stellar output of metals and ionizing radiation 
(Giroux & Shapiro 1996, and references therein). Exam- 
ining the connection between these two effects is critical 
for several problems in cosmology, including the degree to 
which a single population of early stars could have influ- 
enced both reionization and IGM metal pollution, and the 
accurate calculation of the IGM metallicity, Zigm, from 
observations of metal line absorption, which depends on 
the assumptions made for the ionizing photon background 
(e.g., Giroux & Shull 1997). These issues are also per- 
tinent to other fields that are constrained by the stellar 
history of the baryons, such as theories for baryonic dark 
matter candidates (Fields et al. 1998). 

Current spectroscopic data imply that H I reionization 
may occur not far beyond z ~ 6 (Becker et al. 2001), and 
that He II reionization occurs at z ~ 3 (Kriss et al. 2001, 
and references therein). However, recent data suggest a 
higher level of ionizing radiation at z > 6-9 than would 
occur in a scenario with "late" H I reionization at z ~ 6. 
These include observations of: the microwave background 
from the WMAP satellite, indicating a high Thomson op- 
tical depth of - 0.17 ± 0.04 (Kogut et al. 2003) or early 
reionization at z ~ 17 ± 4; the evolution of Ly-limit sys- 
tems and the IGM ionizing background (Miralda-Escude 
2003); and the observed redshift evolution of the Fe/Mg 



ratios in QSO emission-line regions (Hamann & Ferland 

1999) . Although these data are by no means conclusive, 
they indicate that the level of stellar activity at z > 9 
is potentially significant, and that the reionization his- 
tory of the IGM could be relatively complex, with periods 
of extended or double reionization for H I and/or He II 
(Venkatesan et al. 2003; Cen 2003; Wyithe & Loeb 2003; 
Hui & Haiman 2003). Given the current best-fit cosmo- 
logical parameters and the associated small-scale power 
available to the reionizing sources (Spergel et al. 2003), an 
unusually high conversion efficiency of baryons to ionizing 
radiation may be required (Somerville et al. 2003; Ciardi 
et al. 2003). 

At present, it is believed that almost all of the baryons 
at z > 3 reside in the IGM, with an average enrichment 
of Zigm ~ 10" 2 5 Z Q at z ~ 3 (Ellison et al. 2000). This 
is detected primarily through C IV absorption down to 
the lowest column density systems, corresponding to the 
"true IGM"; there is also some evidence for widespread 
O VI enrichment in this IGM component (Schaye et al. 

2000) . From this, one can roughly estimate the fraction 
of baryons that went into the first stars that made the 
IGM metals, and the associated number of ionizing pho- 
tons per baryon (Miralda-Escude & Rees 1997; Haiman & 
Loeb 1997). However, there are two potential problems 
that we face in trying to constrain the effects of high- mass 
stars consistently through the degree of IGM carbon en- 
richment. First, the photons relevant for reionization come 
from the massive stars in the stellar initial mass function 
(IMF; stellar masses > 10 M Q ), while carbon is produced 
dominantly by longer-lived intermediate-mass stars (IMSs; 
~ 2-6 M Q ). These stellar products may not be mutually 
constraining if the IMF was different in the past. Second, 
for burst-driven star formation (not continuous), most of 
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the carbon is produced ~ 0.5 Gyr after the starburst, an 
order of magnitude in time after the massive stars' Type 
II SNc, which occur on timescales of a few tens of Myr and 
which could drive galactic winds. It then becomes unclear 
how the carbon is expelled from individual halos and dis- 
tributed ubiquitously in the IGM by z ~ 3-5 (Songaila 
2001). 

In this Letter, we present calculations of the stellar ion- 
izing efficiency as a function of stellar metallicity, and 
of the generated metals. We specifically examine the re- 
lation between the integrated ionizing radiation and the 
output in carbon and oxygen as: (1) these metals appear 
to be widely detected in the IGM, but are primarily the 
products (and hence probes) of different parts of the IMF, 
and (2) carbon and oxygen are easily observed at typical 
IGM temperatures at z ~ 2-4. The point that carbon is 
primarily the product of IMSs while massive stars dom- 
inate the ionizing radiation is not new to this work, and 
has been demonstrated extensively in the past (see respec- 
tively, e.g., Iben & Truran 1978, and Tumlinson et al. 2003 
and references therein) . Our main goals here are to investi- 
gate the implications of these stellar trends for cosmology, 
particularly for the reionization and metal enrichment of 
the high-z IGM. We focus on stellar nucleosynthetic ar- 
guments, from which we can derive the minimum number 
of ionizing photons that must have been generated in as- 
sociation with the observed levels of IGM enrichment at 
z ~ 2 — 5, including cases with the significant metal yield 
from IMSs. 



2. RELATING STELLAR RADIATION AND METAL YIELDS 

The first generations of stars are expected to be metal- 
free in composition (Pop III), and rely more heavily on 
the p-p chain initially than on the more efficient CNO 
cycle for their thermonuclear fuel source. Consequently, 
they are hotter and emit significantly harder ionizing radi- 
ation relative to their finite- counterparts (Bromm et al. 
2001; Tumlinson et al. 2003; Schaerer 2002, 2003). These 
objects could play an important role in either the single 
or extended reionization of H I and He II. The gain in 
ionizing radiation with metal-free stars depends on the as- 
sumed IMF. Recent theoretical studies indicate that the 
primordial stellar IMF may have been top-heavy (Abel 
et al. 2000; Bromm et al. 2002), leading predominantly to 
stellar masses > 100 M©. 

As noted earlier, for a Salpeter IMF, all of the ion- 
izing radiation and most of the metals, particularly the 
alpha-elements (e.g., oxygen), originate from the short- 
lived massive stars. With decreasing of the stellar pop- 
ulation, the massive stars generate greater ionizing radi- 
ation and lower post-SN metal yield (Heger & Woosley 
2002; Woosley & Weaver 1995). However, the opposite 
trend holds for the IMSs: a star of lower Z+ has greater 
metal yield (van den Hoek & Groenewegen 1997). Such 
constraints on stellar activity therefore alter with the av- 
erage Z+ of starforming regions and hence with redshift. 

As a diagnostic of the relation between the generated 
ionizing photons and metals from stars, Madau & Shull 
(1996) introduced the conversion efficiency, ?7Lyc, of en- 
ergy produced in rest mass of metals, Mzc 2 , to the energy 

4 More recent unpublished calculations indicate that »7Lyc ~ 0.003 for 



released in the H-ionizing continuum, defined as follows: 
VLyc = E(hv > 13.6eV)/(M z c 2 ) (1) 
Madau & Shull (1996) found r)L yc to have a relatively 
constant value of 0.002 4 for solar values of Z+, and to be 
roughly independent of the stellar IMF's slope. This quan- 
tity can be used to calculate the number of ionizing pho- 
tons per baryon in the universe generated in association 
with the observed IGM metallicity, N^ yc /N-t,, defined as 
the ratio of the energy of the Lyc photons per baryon to 
the average energy of a Lyc photon (Miralda-Escude & 
Rees 1997). For solar-Z* stars where r/^ yc = 0.002, we 
have: 

N Lyc /N h = (E Lyc /N h )/(E Lyc ) 

= (VLyc(Zi GM )m p c 2 ) I (E Lyc ) 

~ 0.002 x 10~ 4 x (1 GeV/20eV) - 10. (2) 
This implies that, on average, 10 stellar ionizing photons 
per baryon were generated in association with the syn- 
thesis of the observed Zigm ~ 10 -4 , a reasonable value 
for a number of reasons. Given the observed decline in 
the space density of bright QSOs at z > 3, high- redshift 
star formation is thought to have played a significant, if 
not dominant, role in reionization. Furthermore, although 
A^Lyc/A^b need only equal about unity for reionization, the 
effects of high-z recombinations in the IGM and individual 
halos likely boost the required N-^ yc /Nh to values close to 
about 10 (see Somerville et al. 2003). 

However, 7}h y c is strongly dependent on Z+ and the IMF 
mass range, an effect that has not been studied until re- 
cently (Schaerer 2002). Thus, the above calculation of 
?7Lyc does not always lead to N^ yc /N^ ~ 10. Pop III stars, 
having no metals, are hotter and emit harder radiation. 
When combined with the trend of lower metal yields with 
decreasing Z+, we would expect a net decline in rj-^ yc with 
rising Z*. This has been recently demonstrated for a few 
specific cases by Schaerer (2002). However, the significant 
contribution by IMSs to the overall metal yield, partic- 
ularly through carbon, has been neglected by past works 
(Schaerer 2002; Madau & Shull 1996), which assumed that 
massive stars were the dominant sources of the metals. In- 
cluding the IMS contribution has a non-negligible impact 
on the values of ?7Lyc as well as on N^ yc /N^, as we show 
below. 

A concern in relating the ionizing photon and nucleosyn- 
thetic output of stellar populations is finding an overlap 
amongst the input Z+ and assumptions in the literature's 
calculations. At present, the ionizing photon spectrum for 
metal-free massive stars appears to be consistent amongst 
several works to within ~ 10% (Schaerer 2003; Tumlin- 
son et al. 2003; Bromm et al. 2001). However, the stellar 
yields are not entirely convergent yet, with variations of 
up to factors of a few (see, e.g, Limongi & Chicffi 2002). 
We consider three values of Z+: 0, 0.001 and 0.02 (i.e., 
solar), and two separate mass ranges, 1-100 M© and 100- 
1000 Mq. From here onwards, we refer to the former mass 
range as a present-day IMF and to the latter as very mas- 
sive stars (VMSs) or, equivalently, a top-heavy IMF. Al- 
though such clean distinctions may not exist for the first 
stars, we do this in order to differentiate the signatures of 
stars in these two IMFs. The metal yields are taken from 
van den Hock & Groenewegen (1997) for stellar masses < 

'* = Zq (J. M. Shull 2003, private communication). 
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8 M Q and from Woosley & Weaver (1995) for the 8-40 M 
mass range. Note that for Z* = 0.001 massive stars, we 
use the Z = 0.1Z Q case from Woosley & Weaver (1995). 
Since there are no detailed yields to date for Z* = IMSs, 
we use those for Z* = 0.001 IMSs, which are not sub- 
stantially different (A. Chicffi 2002, private communica- 
tion). For Z* = VMSs, only the stars of mass 140-260 
Mq avoid complete collapse into a black hole (Hcgcr & 
Woosley 2002) and contribute to the nucleosynthetic out- 
put. We note that the term yield in this paper connotes 
the total ejected mass in individual elements or in met- 
als, and not the net yield after accounting for the original 
metal composition of the star. 

We do not consider main-sequence mass loss in any 
of our cases below; this is unlikely to be important for 
low-Z+ stellar populations (see Tumlinson ct al. 2003 on 
this point). We assume the stellar IMF to have the 
form, <j>(M) = 4>oM- a , and the IMF slope, a, to be the 
Salpeter value of 2.35, since small slope variations have 
been shown to make little difference (Giroux & Shapiro 
1996; Madau & Shull 1996; Schacrer 2002, 2003). All the 
cases here are normalized over their respective mass ranges 
as, JdMM<j>(M) = 1. 

For the ionizing spectra, we use Bromm et al. (2001) 
and Schaerer (2003) for a top-heavy IMF, Tumlinson et al. 
(2003) for 1-100 M , Z* = stars, and Leitherer et al. 
(1995) for non-zero Z* cases. We note that the time- 
integrated ionizing photon number from VMSs happens 
to be roughly the same as that from Z* = stars in an 
IMF over 1-100 Mq, despite the greatly boosted ionizing 
rate from VMSs, owing to their brief lifetimes. The stel- 
lar emission rate of ionizing photons drops to less than 
1% of its initial value after 30 Myr for a present-day IMF 
(Venkatesan et al. 2003), but takes only about 3.5 Myr for 
a top-heavy IMF (Schaerer 2003). The cause is directly re- 
lated to the lifetimes of the longest-lived (metal-free) star 
that is relevant for reionization in each of the IMFs: a 
10 Mq star (~ 2 xlO 7 yr) versus a 100 Mq star (~ 3 
xlO 6 yr). The net difference in the time- integrated ioniz- 
ing radiation is a factor of only ~ 1.4 (Bromm ct al. 2001; 
Schaerer 2002, 2003). 

Lastly, for the conversion from 77L yc to Nj^c/N^, we as- 
sume that Zigm ~ 10 _2,5 Z Q , with the correspondingly 
scaled levels of carbon and oxygen with respect to their so- 
lar abundances 5 (Shull 1993, and references therein). We 
take (E Lyc ) to be 21 eV for Z* = Z Q and Z* = 0.001 stars, 
27 eV for Z* = stars in a 1-100 Mq IMF, and 30 eV for 
Z* = stars in a top-heavy IMF (Schaerer 2003). 

3. RESULTS 

We first show that carbon and oxygen are mostly the 
products of IMSs and massive stars respectively, and that 
IMSs arc consequently significant contributors to cosmo- 
logical metal generation. Figure 1 displays the total IMF- 
weighted yields for metals, 12 C and 16 as a function of 
stellar mass for VMSs and for three values of Z* for a 
present-day IMF. Using the total ejected mass, Mj, in each 
element or metals as a function of stellar mass M, the y- 
axis in the figure is calculated by weighting each Mj with 
the IMF as Mj0(M). Using the yields and ionizing spectra 

5 Note that the solar abundance values of some metals have been sigr 
we note that if we use, e.g., Holweger (2001), the values of Nj^ yc /N^ i 
for oxygen. 



as described in §2, we calculate ?7L y c an d Ni, yc /Nb m Table 
1 for several cases: (a) as a function of Z+, (b) as defined 
with respect to total metal, 12 C, and 16 yield, and (c) 
with metal yields included from three mass ranges (1-100 
Mq, 8-100 Mq, and 100-1000 Mq). In the second case, 
the conversion from r/Lyc to A^yc/Vb involves two factors 
that will cause N^ yc /N\, to vary within each column in the 
table: r]^ yc and Ztqm, which will both vary at fixed Z* de- 
pending on whether we are considering rj-^ yc and -/VL yc /iVb 
with respect to total metals, 12 C, or 16 0. The last case is 
intended to evaluate the impact of the metal contribution 
from various regions of the stellar IMF, particularly from 
IMSs, on 7y Lyc and N Lyc /N h . 

The results in Table 1, although subject to the uncer- 
tainties in stellar modelling, reveal several trends of inter- 
est to high-z IGM studies. These may be summarized as 
follows. 

First, ?7Lyc can vary significantly with Z* and with the 
element with respect to which it is defined. As antici- 
pated for the reasons outlined in §2, there is a strong gain 
in ?7Lyc and N^yc/N^ as Z* decreases, with Z = stars in 
a present-day IMF being up to 10-20 times more efficient 
at generating ionizing radiation per unit metal yield than 
the Z = Z© stars in Madau & Shull (1996). 

Second, although the values of ryL yc approximately 
match those of, e.g., Schaerer (2002), when the yields from 
IMSs are excluded, we see that IMSs are in fact a sub- 
stantial source of metals, particularly carbon, and should 
be included in such calculations. One possible exception 
might be at very high redshifts, when IMSs may not have 
had the time yet to eject their nucleosynthetic products. 
In this case, the original definition (Madau & Shull 1996) 
of rjLyc which does not include IMS metal yields would be 
valid, and 7?L yc may be stated to increase with decreasing 
Z* and/or increasing redshift. 

Third, ?7L y c,o is the most sensitive to Z+, increasing 
strongly from solar- Z* to metal-free stars, for a present- 
day IMF. This directly results from the following three 
trends: (1) 16 O and ionizing radiation are produced dom- 
inantly by the IMF's massive stars, (2) with decreasing 
Z*, the production of 16 O decreases while that of ionizing 
photons increases, and therefore, (3) their ratio r?L y c,07 i n 
comparison to ??L y c,z or ryi^c.c, shows the strongest in- 
crease with declining Z* for stars in a present-day IMF. 

Fourth, the value of rihyc directly impacts the fraction of 
baryons required to be converted into a first stars popula- 
tion in order to influence reionization (e.g., Bromm et al. 
2001; Miralda-Escude & Rees 1997). One possible impli- 
cation of the first result above is that a factor of 10-20 
fewer baryons need be part of such a population for them 
to be relevant for reionization alone (which requires only 
Ni,y C /Nb <~ 10), if such stars are not required to account 
for Zigm- In this case, reionization by Pop III stars in a 
present-day IMF occurs late (since they should generate 
fewer ionizing photons than indicated by Table 1), and 
such stars would therefore make a negligible contribution 
to the carbon detected in the high-z IGM, particularly at 
z > 6. An alternate interpretation of Table 1, assuming 
that the IGM is relatively uniformly enriched by volume, 

ficantly revised in recent years. Although this is currently unresolved, 
this work will increase slightly for carbon and decrease by about 30% 
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is that Z = stars (1-100 M ) created the IGM metals 
and consequently, reionization may occur early. This is be- 
cause the values of Ni, yc /Nb (~ 30-220) generated by these 
stars in association with the detected Zjqm well exceeds 
that required for a late H I reionization at z ~ 6. If clus- 
ters of such stars are responsible for both the reionization 
and metal enrichment of the IGM, they should be easily 
detected by future missions such as the JWST (Tumlin- 
son et al. 2003, 2001), and the amount of star formation 
at z > 9 is potentially significant (§1). 

Lastly, despite their high ionizing photon rate, VMSs 
are less efficient at generating total ionizing radiation per 
unit metal mass than are metal-free stars in a present- 
day IMF, owing to the greatly boosted metal yield from 
a top-heavy IMF. (Both of these stellar populations, how- 
ever, have approximately equal ?7Lyc,c> with the difference 
being largest when IMSs are excluded, or equivalently at 
z > 6.) Therefore, VMSs are not necessarily preferred as a 
more efficient source of ionizing radiation at early epochs 
on nucleosynthetic grounds in association with a given de- 
tected level of IGM enrichment, and possibly based on 
reionization requirements as well (see, e.g., Somerville & 
Livio 2003). Furthermore, recent studies (e.g., Schneider 
et al. 2002) indicate that the transition from a top-heavy 
to a present-day IMF occurs at gas metallicities of about 
10 _4 Z Q . If this were true, then we can scale the numbers 
for VMSs in Table 1 down by 1.5 orders of magnitude, 
and derive that VMSs may produce only ~ 0.35 ionizing 
photons per baryon for an IGM metallicity of 10~ 4 Z Q be- 



fore they pollute their environments and cease forming. 
We emphasize, however, that VMSs may still be impor- 
tant for reionization alone, given the right combination of 
conditions, owing to their high ionizing rates. 

The conversion between ?7Lyc and N-^yc/N^ does not ac- 
count for the relative propagation timescales of metals and 
photons from starforming regions to the IGM, which will 
impact the evolution of IGM abundance ratios between in- 
dividual elements, e.g., [C/O]. Furthermore, the definition 
of ?7Lyc implicitly assumes the instantaneous creation of the 
metals. This is not necessarily true at high redshifts, e.g., 
to generate sufficient carbon at z ~ 3-6 requires the star 
formation to have occurred at z > 6-9. As we noted ear- 
lier, this becomes particularly problematic for the trans- 
port of the carbon to the IGM through SN-drivcn winds. 
We explore the cosmological relevance of these timescale 
issues in more detail in a forthcoming paper. 
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Fig. 1. — IMF-wcightcd total ejected masses of metals (solid lines), carbon ( C; dashed lines) and oxygen ( O; dashed-dotted lines) as a 
function of stellar mass for a Salpeter IMF. Top panel: = VMSs in the mass range 100-1000 Mq . Lower three panels: stars of mctallicity 
0, 0.001 and Z e = 0.02 in the mass range 1-100 Mq. 



Table 1 

Ionizing Efficiency as a Function of Stellar Metallicity 





z* = o 


= (1-100 Mq) 


Z* = 0.001 (1-100 Mq) 


Z± = Zq 


(1 100 Mq) 




(10 2 - 10 a M Q ) 


M z (> IM Q ) M z (> 8M ) 


M z (> IMq) M z (> SMq) 


M z {> IMq) 


M z (> 8M Q ) 


»7Lyc,Z (Nhyc/N h ) 
Vhyc,C {N L y C /N h ) 

?7Lyc,0 (N Lyc /N h ) 


0.005 (11) 
0.098 (32) 
0.01 (10) 


0.035 (81) 0.08 (186) 
0.086 (31) 0.48 (173) 
0.15 (172) 0.2 (223) 


0.006 (18) 0.008 (24) 
0.025 (12) 0.088 (41) 
0.013 (18) 0.013 (19) 


0.003 (9) 
0.02 (9) 
0.005 (8) 


0.004 (13) 
0.048 (22) 
0.007 (10) 



Note. — The conversion efficiency of rest mass in metals, carbon and oxygen to ionizing radiation and the associated number of ionizing 
photons per baryon, as a function of stellar metallicity, Z*, and the mass range of the stellar population. Cases are also shown with metal 
contributions from different regions of the stellar IMF. All models assume a Salpeter slope for the IMF. 



